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ABSTRACT 

We present simple analytic approximations for the linear and fully evolved nonlinear mass power spec- 
trum for spatially flat cold dark matter (CDM) cosmological models with quintessence (Q). Quintessence 
is a time evolving, spatially inhomogeneous energy component with negative pressure and an equation of 
state WQ < 0. It clusters gravitationally on large length scales but remains smooth like the cosmological 
constant on small length scales. We show that the clustering scale is determined by the Compton wave- 
length of the Q- field and derive a shape parameter, Pg, to characterize the linear mass power spectrum. 
The growth of linear perturbations as functions of redshift, wq, and matter density f2m is also quantified. 
Calibrating to iV-body simulations, we construct a simple extension of the formula by Ma (1998) that 
closely approximates the nonlinear power spectrum for a range of plausible QCDM models. 

Subject headings: cosmology : theory - dark matter ~ large-scale structure of universe - methods: 
analytical 



1. INTRODUCTION 

Quintessence offers an alternative to the cosmological 
constant (A) as the missing energy in a spatially flat uni- 
verse with a sub-critical matter density Om (Caldwell et 
al. 1998 and references therein). It is an energy compo- 
nent which, similar to A, has negative pressure and there- 
fore a negative wq in the equation of state pg = wq pq . 
However, unlike A for which w — —1, quintessence is 
time evolving and spatially inhomogeneous, and wq can 
have a range of values. The observational imprints of the 
quintessence therefore differ from those of the commonly 
studied ACDM cosmology (e.g., Wang et al. 1999). 

In this Letter, we study spatially flat QCDM models in 
which the cold dark matter and Q-field together make up 
the critical density (i.e. Jim -I- ^q = 1). The quintessence 
is modelled as a scalar field that evolves with a constant 
equation of state wq. It drives the cosmological expansion 
at late times, influencing the rate of growth of structure. 
Fluctuations in Q behave as an ultra-light mass scalar 
field: on very large length scales the quintessence clus- 
ters gravitationally, thereby modifying the level of cosmic 
microwave background temperature anisotropy relative to 
the matter power spectrum amplitude (in addition to a 
late-time integrated Sachs- Wolfe effect); on small length 
scales, fluctuations in Q disperse relativistically and the 
Q-field behaves as a smooth component. 

We investigate the effects of the quintessence on the 
spectrum and time evolution of gravitational clustering in 
both the linear and nonlinear regimes. We propose sim- 
ple, analytical fitting formulas for both the linear and fully 
evolved nonlinear power spectrum of matter density fluc- 
tuations in plausible QCDM models. For the linear power 
spectrum (§ 2), we introduce a simple parameter Pq de- 



rived from the Compton wavelength of the Q-field to char- 
acterize its shape. This parameter determines the length 
scale above which the Q-field can cluster gravitationally, 
and is reminiscent of P^ derived from the free stream- 
ing distance of hot neutrinos in cold+hot dark matter 
(C-I-HDM) models by Ma (1996). For the nonlinear power 
spectrum (§ 3), we examine the validity of the simple lin- 
ear to nonlinear mapping technique that has been success- 
fully developed for scale-free, CDM, C+HDM, and ACDM 
models (Hamilton et al. 1991; Jain, Mo, & White 1995; 
Peacock & Dodds 1996; Ma 1998). We present a simple ex- 
tension of the analytical formula of Ma (1998) that closely 
approximates the QCDM nonlinear power spectrum com- 
puted from a set of A^-body simulations. 

The formulas presented in this Letter are essential for 
gaining physical insight into the effects of the quintessence 
on gravitational collapse and for performing rapid predic- 
tions of observable quantities in the linear as well as non- 
linear regimes in plausible QCDM models. 

2. LINEAR POWER SPECTRUM 

We use the conventional form to express the linear 
power spectrum^ for the matter density perturbation 5m 
in QCDM models: 

P(fc,a)=AQfc"r|(fc) f^ , (1) 

\9QflJ 

where Aq is a normalization, k is the wavenumber, n is 
the spectral index of the primordial adiabatic density per- 
turbations, and Tq is the transfer function which encapsu- 
lates modifications to the primordial power-law spectrum. 
The function qq is the linear growth suppression factor, 



Defined so that the two-point correlation function is ^(r) = 47r / k dk P{k)sin{kr)/(kr). 



qq ~ D/a, where D is the standard hnear growth fac- 
tor for the matter density field in QCDM models, and 
9Q,o = 5Q(a = 1) denotes its value at the present day with 
scale factor a — 1. We discuss each piece of equation (|l]) 
in turn. 
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FIG. 1 — Ratio of the transfer functions, Tqa = Tq/Ta, at 
the present day for six pairs of flat QCDM and ACDM models. 
The solid (for f2ni = 0.4) and dashed (fim = 0.6) curves are 
computed from the Boltzmann integrations and illustrate the 
dependence on the matter density parameter flm. For a given 
Qm, the three curves illustrate the dependence on the equation 
of state: wq = —1/3, —1/2, and —2/3 from top down. The 
dotted curves show the analytic approximation given by equa- 
tion (H)-d4|). Note that Tqa deviates from unity only on very 
large length scale {k <^ 0.01 h Mpc~^) above which the Q-field 
can cluster spatially. 



First we examine the transfer function for the matter 
density field. To isolate the effects of quintessence, we 
find it convenient and illuminating to compare a pair of 
QCDM and ACDM models that have the same set of cos- 
mological parameters and differ only in wq (recall w ~ —1 
for ACDM). We define the relative transfer function for a 
pair of such models to be Tqa = Tq/Ta- For Ta, we follow 
the convention and set the arbitrary amplitude of Ta to 
unity as fc ^ 0. The form of Ta is well known and various 
fitting formulas have been published (e.g., Bardeen et al. 
1986; Efstathiou, Bond, & White 1992; Sugiyama 1995). 
More complicated fits with higher accuracy have also been 
developed for higher baryon ratios {flh/^m <; 20%) and 
for the features due to baryonic oscillations and damping 
(e.g., Bunn & White 1997; Eisenstein & Hu 1998). 

The transfer function Tq for QCDM models resembles 
Ta for ACDM but with one key difference. The linear 
matter density field, Sm, evolves according to the equa- 



tion 6m + 'iHSm = 47rG(pm(5m + SpQ + SSpg), where 
H = d/a and the dots denote differentiation with re- 
spect to proper time. On small length scales, the Q- 
field is smooth {Spg, 6pQ ^ Spm) and we recover the fa- 
miliar equation for the evolution of Sm (Caldwell et al. 
1998). On very large length scales, however, the Q-field 
clusters and contributes to the energy density and pres- 
sure perturbations. The result is a different growth rate 
for 6m on large and small scales once the quintessence 
starts to dominate the cosmological energy density. We 
can determine the characteristic scale separating these 
two regimes by examining the linear equation for the Q- 
field: SQ + 3H6Q + [P + Vqq)6Q = (5„[(1 + wq)pqYI^ 
(where V is the Q-field potential, V,qq = (PV/dQ^, and 

Pq = Q^/2 + V). We see that 6Q itself behaves as a 
scalar field with an effective mass (V^qq)^'^ and a Comp- 
ton wavenumber of kg ~ (Kqq)""^^^- On small length 
scales (i.e. k ^ kg), the amplitude of 6Q and hence 6pQ 
is damped and does not enter the evolution equation for 
6m- On large scales {k <^ kg) 6Q grows, so the Q-field 
clusters and in turn affects the evolution of 6m - 

The change in the behavior of 6m near k '^ (V.Qqy''^ as 
a result of differing Q-clustering properties is illustrated 
in Figure 1 for Tqa vs. k for a range of wq and f2i„. 
We have chosen to normalize Tqa to unity at the high-fc 
end because both the Q-ficld and the cosmological con- 
stant are spatially smooth on these scales. The clustering 
property of Q is reminiscent of the case of massive neutri- 
nos in C-I-HDM models, which cannot cluster appreciably 
below the neutrino free-streaming scale but can cluster 
with the same amplitude as the cold dark matter on large 
scales. Analogous to the shape parameter F^ that was 
introduced to model the neutrino streaming distances in 
C-I-HDM models (Ma 1996), we introduce a new shape 
parameter Tq here to characterize the feature in Figure 1 
in QCDM models. For a constant equation of state, wq, 
one can show that Vqq — 67rG(l — wq){2p + p + wgp) , 
where p and p are the total energy density and pressure. 
We approximate 



kQ = TQh = 2y%QQ 
3H 
c 



'J i_i I 

^Jil~WQ)[2 + 2wQ-WQnn,ia)] (2) 



and we use a simple ratio of polynomials to express the 
relative transfer function: 



TQA{k, a) = 7— = -— J , 

Ta 1 + aq^ 



k 



q = 



^Qh 



(3) 



where k is in Mpc^^, and a is a scale- independent but 
time-dependent coefficient that quantifies the relative am- 
plitude of the matter density field 6m on large and small 
length scales. We find a well approximated by 

a^i-wgY, (4) 

s= (0.012 - 0.036 WQ - 0.017/u;q)[1 - n^{a)] 
+ (0.098 -t- 0.029 WQ - 0.085/wq) In an(a) , 

where the matter density parameter is flmio-) — f^m/[^in + 
(1— rim) a~'^™'3] ^ which reaches the value flm at the present 
day a — 1- Figure 1 illustrates the close agreement (with 



errors < 10%) between the approximations given by equa- 
tions (y) - (y) and the exact results from numerical inte- 
grations of the Boltzmann equations. 
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FIG. 2 — Ratio of the growth suppression factors, qqa = 
qq/qa, as a function of the scale factor a (top) and the equation 
of state WQ (bottom; at o = 1) for various pairs of flat QCDM 
and ACDM models. The solid and dashed curves in both pan- 
els are for f7m = 0.4 and 0.6, respectively. In the top panel, 
each set of curves corresponds to wq — —2/3, —1/2, —1/3, and 
— 1/6 from top down. 



Next we examine the linear growth suppression factor of 
the density field in equation (|l|) . This function is well stud- 
ied for ACDM models (Heath 1977; Lahav et al. 1991). 
An empirical fit is given by g\ = 2.5riin(a){^m(a)^'^ — 
1 + an(a) + [1 + f^m(a)/2][l + (1 - f^m(a))/70]}-i and is 
accurate to ~ 2% for 0.1 < flm < 1 (Carroll, Press, & 
Turner 1992). This formula unfortunately cannot be gen- 
eralized to QCDM models by simply replacing (1 — Sim) — > 
(1 — S7in)a^'^(^+"'"3^. Instead, we propose the following for- 
mula to approximate the ratio of the QCDM and ACDM 
growth factors: 



5QA = 



_9Q_ 
5A 



~^qY 



(5) 



i=-(0.255 + 0.305 wq + 0.0027/uiq)[1 - n„,{a)] 
- (0.366 -t- 0.266 wq - 0.07 /wq) In n„,{a) 

accurate to 2% for 0.2 <; 17^ < 1 and —1 < wq ^ —0.2. 
Figure 2 illustrates the dependence of ggA on time, wq, 
and Qrn- The growth is evidently slower in models with 
less negative wq for fixed fim- This is because the en- 



ergy density in the Q-field dominates over that in matter 
at an increasingly earlier tinre as wq is varied from —1 
to 0; the growth of gravitational collapse therefore ceases 
earlier and results in a smaller value for gqA- (It is some- 
times useful to study the instantaneous growth rate of (5,„ , 
/ = d\ogSm/d\og a. See Wang & Steinhardt (1998) for a 
fitting formula for /.) 

The remaining component in equation dlf) to be speci- 
fied is the normalization Aq. It can be chosen by fixing 
the value of ag, the rms linear mass fluctuation within 
a top hat of radius 8h~^ Mpc, or by fixing to COBE 
results. For the latter we follow Bunn & White (1997). 
In the case that the temperature anisotropy is due to 
primordial adiabatic density perturbations with spectral 
index n, we write Aq = 5'jj{c/ H^))"'^'^ / {An) , where 5h 

5^,-1 



2 X lO-^ao"' {a 



)ci+C2 lnO„, gj^p [^gj-^ -l)+ Ci{n - 1)2] , 

ofQ — a{a — 1) of equation (0), ci = 
-0.789|wq|O-o^5'*-''-2ii1"^"qI, c2 = -0.118 - 0.0727u;q , 
C3 = -1.037, and ca = -0.138, for -I <, wq <, -0.2. In 
the case of tensor perturbations, the primordial amplitudes 
follow the inflationary relation At = 8(1 — n)As ■ The ra- 
tio of £ = 10 multipole moments, rio = Cj^/Cfg, is given 
by no « 0.48(1 - n)[l + 0.1(1 - n)(8 + Iwq)^^'' + 3](1 - 
r2Q/x)si«(^«/"') , where 510(2/) = 0.18 -{- 0.84y2 , and 
X = 0.75[1 - 0.66u;q -I- 1.66w| - 0.5(1 -f wqf] . One 
can rescale Aq by Aq — *■ Aq/{1 + rio) to accommodate 
the effect of tensors on the normalization. 

3. NON-LINEAR MASS POWER SPECTRUM 

In this section we examine if the simple linear to non- 
linear mapping technique initiated by Hamilton et al. 
(1991) can be extended to QCDM models. The basic ap- 
proach is to search for a simple expression for the func- 
tion Ani(fc) = /[Ai(fco)] that relates the linear and nonlin- 
ear density variance A(fc) = ATTk^P{k). Note that Ani 
and Al are evaluated at different wavenumbers, where 
fco — fc(l -I- Ani)~^/^ corresponds to the precoUapsed scale 
oik. The strategy is to combine analytical clustering prop- 
erties in asymptotic regimes with fits to numerical simu- 
lation results. This recipe has been successfully developed 
for scale-free models with a power-law P{k) (Hamilton et 
al. 1991; Jain et al. 1995), flat CDM and ACDM models 
(Jain et al. 1995; Peacock & Dodds 1996, PD96 hereafter; 
Ma 1998), and flat C-I-HDM models with massive neutri- 
nos (Ma 1998, Ma98 hereafter). 

We investigate if the PD96 and Ma98 formulas pro- 
posed for ACDM models can be easily extended to QCDM 
models. These two formulas incorporate the time depen- 
dence of the mapping in different ways, but they share 
the feature that the dependence on parameters f2m and 
r^A enters only through the linear growth factor g. In 
order to test the application of this method to QCDM 
models, we have performed A'^-body simulations for three 
values of wq: -2/3,-1/2, and —1/3, each with several 
different realizations. These three values should be suffi- 
cient since extensive tests of wq — —1 (i.e. ACDM mod- 
els) have already been carried out in PD96 and Ma98. 
We restrict our attention to wq < —1/3 and cosmolog- 
ical parameter ranges that are in concordance with ob- 
servations (Wang & Steinhardt 1998; Wang et al. 1999). 
Specifically, (f^^, fiq, fib, h) = (0.4, 0.6, 0.047, 0.65) for the 
wq — —2/3 and —1/2 models, and (fini, ^q, ^b, ^) = 



(0.45,0.55,0.047,0.65) for the wq = -1/3 model. The TV- 
body code used is a parallel version of the particle-particle 
particle-mesh algorithm (Bertschinger & Gelb 1991; Fer- 
rell & Bertschinger 1994). Each simulation uses 128'^ par- 
ticles in a box of comoving volume 100'^ Mpc"^. The Plum- 
mer softening length is 50 kpc comoving, which allows us 
to compute the nonlinear power spectrum in highly clus- 
tered regions with k <^10h Mpc~^ and Ani <; 1000. Since 
the Q-field clusters only on scales much above the box 
size, the presence of the quintessence only affects the ini- 
tial conditions and the evolution of the scale factor a. 
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FIG. 3 — Lower left and upper. The linear and fully evolved 
nonlinear power spectra for the matter density field in QCDM 
models with different equations of state wq (see text for other 
model parameters). In each panel, five redshifts, ^ = 0, 1, 2, 
3, and 4, are shown (top down). The curves are computed 
from: A'^-body simulations directly (thick solid), nonlinear ap- 
proximation by Ma98 (dashed; eq. (6)) and PD96 (dotted), 
and linear theory (thin solid). Lower right: The effective spec- 
tral index, ricff, as a function of as for two QCDM models. 
The dotted line represents a power law and demonstrates that 
dln(neff -|- 3) /din as oc /3 is an excellent approximation. 



Figure 3 compares the linear power spectrum and the 
fully evolved spectrum from both the iV-body runs and 
the approximations of PD96 and Ma98. Five redshifts are 
shown for each of three QCDM models. Overall, we find 
that the PD96 formula works well at z = when the factor 
g in their fornrula is set to g — gg, which is the appropri- 
ate growth factor for the density field for QCDM models 
[eq. (R)]. At earlier times, however, the PD96 formula un- 
derestimates Ani at /c <; 1 h Mpc^^ in the wq — —2/3 



and —1/2 models by up to 30%. We have attempted less 
physically- motived combinations of growth factors (e.g., 
g — a gq and ^a) but did not find a way to make PD96 
fit. We find that the Ma98 formula, 



Anl(fc) 

Ai(fco) 



G 



A, 



3/2 /3 

,9o cTs, 



G{x) = [l + ln(l + 0.5a;)] 



l + 0.02a;4 + cia;Vg3 



1 + C2 a;' 



7.5 



(6) 



can be easily extended to QCDM models. Specifically, we 
propose to keep ci = 1.08 x 10^** and C2 — 2.10 x 10~^ 
used for ACDM in Ma98, but adopt g — gg, which 
is the appropriate QCDM growth factor [eq. (||)], and 

50 - |w;q|^-3|'"«I-°-^^5q,o, where 5q,o = ggia =1). As 
described in Ma98, the parameter /3 in equation (|g) is in- 
troduced to approximate the power-law dependence of the 
effective spectral index rics + 3 in previous work on erg: 
d\n{ncff + 3)/d\nag, oc (3 . We find (3 = 0.83 an excellent 
approximation for all three QCDM models that we tested 
(see the lower- right panel of Figure 3). Other panels of 
Figure 3 illustrate the close agreement (rms errors ^ 10%) 
between TV-body results and equation (H). 



4. SUMMARY 

We have presented simple formulas to approximate both 
the linear and nonlinear power spectra for matter density 
perturbations in viable quintessence cosmological models 
with an equation of state —1 < wq ^ —1/3. Equa- 
tions (p|),(H), and (0) together specify the ratio of the linear 
transfer functions Tq and T\ for the matter density field 
for a given pair of QCDM and ACDM models with the 
same cosmological parameters. Equation (o) specifies the 
ratio of the linear growth suppression factors gq and ^a 
in QCDM and ACDM models. Equation (ph approximates 
the nonlinear mass power spectrum. 

A key difference between gravitational clustering in 
QCDM and ACDM models is that A is spatially smooth 
on all length scales, whereas the Q-field can cluster above 
a certain length scale. We characterize this length scale by 
the shape parameter Tq of equation (0), which is derived 
from the Compton wavelength of the Q-field. The QCDM 
matter power spectrum therefore changes shape at two 
characteristic scales: Tq, and the familiar F oc fluih that 
corresponds to the cross-over from radiation- to matter- 
dominated era. For the QCDM models studied in this 
Letter (i.e. constant wq), the Compton wavelength of the 
Q-field is very large: fcg ~ 0.001 to 0.01 /i Mpc'^ On 
scales of galaxy clusters and below, therefore, the linear 
QCDM power spectrum has identical shape as in the cor- 
responding ACDM model and differs only in the overall 
amplitude by a factor of [Aq / Ai^)[gQf^l gQf^fiY' . This re- 
alization should simplify comparisons between QCDM and 
ACDM models. 

For the fully evolved nonlinear power spectrum, we find 
that PD96 works well at z = but underestimates its am- 
plitude by up to ~ 30% at earlier times. The formula of 
Ma98, on the other hand, can be easily extended to ap- 
proximate the QCDM nonlinear P{k) (with errors < 10%) 
for Wq < —1/3 and redshift up to z w 4. Equation (|^) 
summarizes this result. 
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